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Abstract 

We extract a relatively precise value for the decay constant of the meson by measuring 
B{D+ fi+u) = (4.40 ± 0.66l:g:?|) x 10"^ using 281 pb^^ of data taken on the ^(3770) resonance 

with the CLEO-c detector. We find = (222.6 ± 16.7l|^) MeV, and compare with current 
theoretical calculations. We also set a 90% confidence upper limit on B{D~^ e'^v) < 2.4 x 10~^ 
which constrains new physics models. 
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Purely leptonic pseudoscalar meson decays proceed via the annihilation of the constituent 
quarks into a virtual W via the axial- vector current. The decay rate is proportional to the 
square of the decay constant /, a single number which encapsulates strong interaction dy- 
namics in the decay. The lack of hadrons in the final state allows for precision tests of strong 
interaction theories. Knowledge of decay constants is critical for extracting fundamental 
information about CKM matrix elements. For example, Jb is needed to use measurements 
of BB mixing. Currently, it is not possible to determine fs experimentally, so theoretical 
calculations must be used. The most promising of these calculations involves lattice QCD P], 
though there are other methods. Measurements of pseudoscalar decay constants from charm 
meson decays provide checks on these calculations and can discriminate among different 
models. 

The decay —>■ i^u proceeds by the c and d quarks annihilating into a virtual W^, 
with a decay width |2|] 

T{D^^t.) = ^fl,m^,M^,l^l-^J\Vj' , (1) 

where is the mass, me is the mass of the final state lepton, \Vcd\ is a CKM matrix 
element that we assume to be equal to \ Vus\, and Gp is the Fermi coupling constant. Because 
of helicity suppression, the rate is a function of mj; consequently, the electron mode 
e^v has a very small rate in the standard model. The relative widths are 2.65 : 1 : 2.3 x 10^^ 
for the r"*"//, ii^y, and e'^u final states, respectively. 

The CLEO-c detector is equipped to measure the momenta and directions of charged 
particles, identify them using specific ionization (dE/dx) and Cherenkov light (RICH), detect 
photons and determine their directions and energies j3|. 

In this study we use 281 pb~^ of data produced in e^e~ collisions using the Cornell 
Electron Storage Ring (CESR) and recorded at the tp" resonance (3.770 GeV). This work 
contains our previous sample as a subset and supersedes our initial effort j^. At this energy, 
the events consist mostly of D^D^, D^D , three-flavor continuum qq [q = u,d,s), t^t^, 
and 'yip' events. 

Our analysis strategy is to fully reconstruct the D~ meson in one of six decay modes listed 
in Table H] and search for a —>■ fi'^u decay in the rest of the event. Charge conjugate 
modes are implicitly included throughout the paper. Track selection, particle identification 
(PID), 7r°, Ks, and muon selection cuts are identical to those used in Ref. 0]. We first 
evaluate AE, the difference in the energy of the decay products with the beam energy. The 
AE distributions in all modes are well described by either a Gaussian or the sum of two 
Gaussians, with root mean square (r.m.s.) widths varying from 7 MeV for K^K^tx^ 
to 14 MeV for K+'r-ti-'k^. We select candidates by requiring |A£^| <0.012 - 0.024 

GeV, where the cut in each mode is approximately 2.5 times the r.m.s. width. 

For the selected events we construct the beam-constrained mass, m^c = 
\/-^beam ~ i^iPiYi wherc i runs over the final state particles from the candidate D~ de- 
cay. The resolution in itibc of 2.2-2.4 MeV is better than merely calculating the invariant 
mass of the decay products, since the CESR beam has a small energy spread. The niBc 
distribution for the sum of all D~ tagging modes is shown in Fig. ^ The numbers of tags 
in each mode are determined from fits of the itlbc distributions to a signal function plus a 
background shape. For the latter we use an expression analogous to one first used by the 
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ARGUS collaboration to approximate the correct threshold behavior [5]. For the signal we 
use an asymmetric lineshape because of the tail towards high mass caused by initial state 
radiation Table H] gives the numbers of signal and background events for each mode 
within the signal region, defined as mo — 2.5 CTmBc < "^sc < i^d + 2.0 (TmBc^ where (TmBc 
the r.m.s. width of the lower side of the distribution. 



Mode 




Signal 


Background 




TT 


77387 ± 281 


1868 


K+ir- 




24850 ± 214 


12825 






11162 ± 136 


514 






18176 ± 255 


8976 






20244 ± 170 


5223 


K+K 




6535 ± 95 


1271 


Sum 




158354 ± 496 


30677 



TABLE I: Tagging modes and numbers of signal and background events. 
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FIG. 1: Beam-constrained mass for the sum of fully reconstructed D~ decay candidates. The solid 
curve shows the fit to the sum of signal and background functions, while the dashed curve indicates 
the background. 

Using our sample of D~ candidates we search for events with a single additional charged 
track presumed to be a /i"*". The track must make an angle >35.9° with respect to the beam- 
line, deposit less than 300 MeV of energy in the calorimeter, characteristic of a minimum 
ionizing particle, and not be identified as a kaon. We then calculate 

= (Ebeam ' E^+f - {-p^- - p^+f , (2) 

where Pj-,- is the three-momentum of the fully reconstructed D~. Real —>■ /i+z/ events 
will congregate near zero MM^. 

The MM^ from Monte Carlo simulation has a resolution (two-Gaussian a) of 
0.0235±0.0004 GeV^ consistent among all the tagging decay modes. We check our simu- 
lations by using the D"^ — > KgT^^ decay. Here we choose events with the same requirements 
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as used to search for /i+z^ but require one additional found Ks- The resolution is measured 
to be 0.0233±0.0009 GeV^, consistent with the Monte Carlo estimate. 

In order to restrict the sample to candidate ^'^u events, we impose restrictions on tracks 
and neutral energy clusters in addition to those used to reconstruct the tagging D~. We select 
events with only one additional charged track; events with extra tracks originating within 0.5 
m (radially) of the event vertex are rejected. In addition, we eliminate events having at least 
one extra neutral energy cluster of more than 250 MeV. These cuts are highly effective in 
reducing backgrounds especially from — > tt+tt^ decays, but they introduce an inefficiency 
because the decay products of the tagging D~ can interact in the detector material leaving 
spurious tracks or clusters. To evaluate our cut efficiencies, we use an essentially background- 
free sample of fully reconstructed D~^D~ events. (The method is different here than in our 
original publication, though the results are consistent.) 

To first order, the fully reconstructed D~D^ —>■ X+tt^tt^, K^n^n'^ events can be consid- 
ered the superposition of two D~ K^n^n^, fi^u events. Our procedure is to evaluate 
the cut efficiency in our sample of 1435 events and take the square- root. This gives us the 
efficiency for the D~ K^n^n^ tag sample. We then combine the K^tt^tt" with each 
of the other tags in turn. This method ensures that the number of interactions of particles 
with material and their resulting effects is the same as in the tag sample used for the fi'^u 
analysis. In the sample of fully reconstructed D^D^ events there are no events with extra 
tracks originating within 0.5 m of the main event vertex. The efficiency for rejecting events 
with extra clusters above 250 MeV, averaging over all our tag modes, is (96.1±0.3±0.4)% 
The systematic error arises only because we have analyzed a situation corresponding to two 
overlapping tags rather than one tag plus a muon. In Monte-Carlo simulation the efficiency 
difference is 0.4%, which we assign as a systematic error. 

The MM^ distribution is shown in Fig. |21 We see a peak near zero containing 50 events 
within the interval —0.050 GeV^ to +0.050 GeV^, approximately ±2cr wide. The peak is 
mostly due to — > fi^iy signal. The large peak centered near 0.25 GeV^ is from the decay 
k'^tt^ that is far from our signal region and is expected, since many escape our 
detector. 

There are several potential background sources; these include other modes, misidenti- 
fied D^D events, and continuum including e^e~ jip'. Hadronic sources need to be consid- 
ered because the requirement of the muon depositing less than 300 MeV in the calorimeter, 
while about 99% efficient on muons, rejects only about 40% of pions or kaons as determined 
from a pure sample of K~7t~^ decays. 

There are a few specific decay modes that contribute unwanted events in the signal 
region. Residual 7r"'"7r'^ background is determined from a simulation that uses a branching 
fraction of (0.13±0.02)% and yields 1.40±0.18±0.22 events; the ffist error is due to Monte 
Carlo statistics, and the second is systematic, due mostly to the branching ratio uncertainty. 
We find background from — » T~^iy only when r+ tt^u. Since the T~^iy branching ratio is 
known to be 2.65 times the fi'^u rate from Eq.^ our simulation gives 1.08±0.15±0.16 events, 
where the systematic error arises from our final uncertainty on the fi'^u decay rate. The K tt"'" 
mode (branching ratio of (2.77±0.18)% [^) gives a large peak in the MM^ spectrum near 
0.25 GeV^. While far from our signal region, the tail of the distribution can contribute. We 
measure this background rate directly using D^D^ events. Here we select a sample of single 
tags, either K~7r~^'n'~^7r~ , K~'k^tt^ or K~tt^ and look for events with only two additional 
oppositely signed tracks, one identified as a kaon and one as a pion using the RICH. We 
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FIG. 2: MM^ using tags and one additional opposite sign charged track and no extra energetic 
clusters (see text). The insert shows the signal region for — n'^v enlarged; the defined signal 
region is shown between the two arrows. 



then compute the MM^ ignoring the kaon. The MM^ distribution shows a narrow peak 
near 0.25 GeV^ and three events in the signal region, corresponding to a background of 
0.33±0.19±0.02 events, the systematic error being due to the branching ratio uncertainty. 
(A simulation gives a consistent estimate of 0.44±0.22 events.) 

We have also checked the possibility of other D^D^ decay modes producing background 
with an equivalent 1.7 fb~^ Monte Carlo sample. We evaluate D'^D and continuum back- 
grounds by analyzing Monte Carlo samples corresponding to 0.54 fb^^. To normalize our 
Monte Carlo events to our data sample, we use cr^o;^" ~ "^"^ cTcontinuum = 14.5 nb (9i]. 

No additional background events are found in any of these samples. 

Our total background is 2.81±0.30±0.27 events. The backgrounds from other D^, D^, 
and continuum sources are limited to less than 0.4, 0.4, and 1.2 events at 90% confidence level 
(C.L.), respectively. To account for possible backgrounds from these sources, we add them 
as 32% C.L. (la) values in quadrature for a positive error and therefore add an additional 
^0® event systematic error. 

We have 47.2±7.llo;8 /i+z/ signal events after subtracting background. The detection 
efficiency for the single muon of 69.4% includes the selection on MM^ within ±2a limits, 
the tracking, the particle identification, probability of the crystal energy being less than 300 
MeV, and corrections for final state radiation [ifli]. It does not include the 96.1% efficiency 
of not having another unmatched cluster in the event with energy greater than 250 MeV. We 
also need to account for the fact that it is easier to find tags in /i+i^ events than in generic 
decays by a small amount, (1.5±0.4±0.5)%, as determined by Monte Carlo simulation. 

Our result for the branching fraction, using the tag sum in Table 12 is 



(4.40 ± 0.66tn^) X 10 



(3) 



The systematic errors on the branching ratio are listed in Table |n] (The systematic error 
on the tag sum is estimated from varying the signal and background functions.) 
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Systematic errors (%) 



MC statistics ±0.4 
Track finding ±0.7 
PID cut ±1.0 
MM2 width ±1.0 
Minimum ionization cut ±1.0 
Number of tags ±0.6 
Extra sliowers cut ±0.5 

Background ±0.6, —1.7 

Total ±2.1, -2.5 



TABLE II: Systematic errors on the D"*" ii^u branching ratio. 

The decay constant fo* is then obtained from Eq. using 1.040±0.007 ps as the D'^ 
hfetime 0, and \Vcd\ = 0.2238±0.0029 yT|]. (We add these two small additional sources of 
uncertainty into the systematic error.) Our final result is 

/d+ = (222.6 ± l&.ltll) MeV . (4) 

We use the same tag sample to search for e^i^e- We identify the electron using a 

match between the momentum measurement in the tracking system and the energy deposited 
in the Csl calorimeter as well as insuring that the shape of the energy distribution among the 
crystals is consistent with that expected for an electromagnetic shower. Other cuts remain 
the same. We do not find any candidates, yielding a 90% C.L. limit of B{D^ e^i'e) < 
2.4 X 10~^ , including systematic errors. 

Our measurement of /^n- is much more precise than previous observations or limits |3.E^. 
The theoretical predictions listed in Table UTTl were made prior to this result. The first entry 
is the result from the Fermilab-MILC-HPQCD collaboration that is done with all three light 
quark flavors unquenched, hence nj=2±l [13]. It is about 10% smaller than our result, albeit 
within error. 

The models generally predict fj-)+ to be 10-25% larger than fD+ which is consistent with 
a previous CLEO measurement jl^]. Some non-standard models predict significant rates for 
the helicity suppressed decay — > e^i^ ji^. Our upper limit restricts these models. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent lumi- 
nosity and running conditions. This work was supported by the National Science Foundation 
and the U.S. Department of Energy. 
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